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MEM study of true flattening of free energy and the 9 term*^ * 
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We study the sign problem in lattice field theory with a 9 term, which reveals as flattening phenomenon of the 
free energy density f(6). We report the result of the MEM analysis, where such mock data are used that 'true' 
flattening of f(8) occurs . This is regarded as a simple model for studying whether the MEM could correctly 
detect non trivial phase structure in 9 space. We discuss how the MEM distinguishes fictitious and true flattening. 



1. INTRODUCTION 

Lattice field theory with the 9 term suffers from 
the sign problem. A conventional technique to 
circumvent it is to perform the Fourier transform 
of the topological charge distribution P(Q) in or- 
der to calculate the partition function Z{9) pQ. 
The distribution P(Q) is calculated with real pos- 
itive Boltzmann weight. However, this still causes 
flattening of the free energy, which misleads to a 
fictitious signal of a first order phase transition. 
To overcome this problem requires exponentially 
increasing statistics with lattice volume. 

We consider the issue of flattening in terms of 
the maximum entropy method (MEM) 2 3. In 
our previous paper 0] (referred to as (I) here- 
after), we used mock data constructed by the 
Gaussian P(Q) whose corresponding free energy 
is analytically known to have no singular struc- 
ture at values of 8 for < 9 < it and tested 
whether the MEM would be effective in the case 
with flattening as well as without flattening. In 
the case without flattening, the results of the 
MEM agree with the exact results. In the case 
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with flattening, the MEM gives smoother f{ff) 
than that of the Fourier transform. 

In the present work, we consider an opposite 
case to the above, i.e., consider a model that 
causes real flattening, which will be explained 
later. This mimics a case where non trivial phase 
structure occurs at finite value of (9 < ir) . We 
generate mock data based on the model and ap- 
ply the MEM to them. Our aim is to study (i) 
whether flattening would be properly reproduced 
in this model and (ii) how true and fictitious flat- 
tening are distinguished in the MEM analysis. 

2. MEM AND 'TRUE' FLATTENING 

The partition function Z(9) can be obtained by 
Fourier-transforming the topological charge dis- 
tribution P(Q): 



Z(9) 
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} P(Q), (i) 



where S is the action, and Q{z, z) is the topologi- 
cal charge described by lattice fields z and z. The 
distribution P(Q) obtained from MC simulations 
can be decomposed into two parts, a true value 
and an error of P(Q). When the error at Q = 
dominates, the free energy density could develop 
a flat region in the large 9 region, which misleads 
one into identifying flattening as a signal of a first 
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order phase transition at a finite value of 9 |5l(:i| . 
Instead of dealing with Eq. (JIJ, we consider 
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The MEM is based on Bayes' theorem. In our 
formulation, in order to calculate Z(0), we max- 
imize the conditional probability, called the pos- 
terior probability 



prob(Z(<9)|P(<3),/,a,m) oc e" 



(3) 



where the probability is determined by \ 2 an d 
the Shannon- Jaynes entropy S which includes the 
default model m(9) 



S = 



(10 



Z{6) - m(6) - Z{6) log 



Z{6) 



m{9) 



(4) 



Information / plays the important role that en- 
forces Z{9) to be positive. The parameter a de- 
termines the relative weights of S. The best im- 
age of Z(9) for a fixed a is then given by 
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This is followed by the successive procedures 
in order to obtain a-independent final image 
Z{9) 2 ,'. : (i) Averaging Z^: 

Z{9) = [ dapiob{a\P(Q),I,m)Z < - a \6). (6) 



and (ii) error estimation. 

In (I), we used the Gaussian P{Q) for the MEM 
analysis, which is parameterized as 



P G (Q) ocexp[--Q 2 ] 



(7) 



where c is a constant (=7.42) and V is regarded 
as the lattice volume. In order to address a ques- 
tion whether flattening is due to the statistical er- 
ror or due to the characteristic of the data them- 
selves, and a question whether these two are dis- 
tinguished in terms of the MEM, we consider a 
simple model. Suppose in some lattice theory 
that calculated P(Q) were slightly deviated from 
the Gaussian one only at Q = 0: 



P(Q) = P G (Q)+Pox6 Qfi , 



(8) 
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Figure 1. Partition function calculated by nu- 
merical Fourier transform for V = 50. Data for 
P{Q) are based on Eq. © with p = 4.5 x 10~ 2 , 
4.5 x 1(T 3 , 2.0 x 10~ 3 and 0. 



where po is a constant. The partition function is 
analytically obtained by use of the Poisson sum 
formula as 



Z ps (9) cx 



exp 



-27m) 



-po-(9) 



Because the first term in Eq. © is monotonically 
decreasing function of 9 in the region < 8 < n, 
the second one gives a flat distribution at large 
values of 9. We generate mock data based on 
P(Q) in Eq. JHJ by adding Gaussian noise with 
the variance 6 x P(Q) for each value of Q. 

3. RESULTS 

Figure displays the behavior of Z(9) cal- 
culated by the numerical Fourier transform of 
Eq. JHJ for V = 50 and various values of po; 
p a = 4.5 x 10~ 2 , 4.5 x 10~ 3 and 2.0 x 10~ 3 . The 
parameter S is chosen to be 1/400. Figure ^also 
displays Z(9) for V — 50 and po = 0, which was 
studied in (I) in detail. Clear flattening is ob- 
served in all the cases. When noise is added, the 
partition function for Eq. is given by 
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Figure 2. Z^ a \6) for various values of a. Here, 
V = 50 andp = 4.5 x 1CT 2 . 
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Figure 3. Z(9) for V = 50 and p = 4.5 x 10" 



where B is a normalization constant. When the 
constant term dominates in the large 9 region, 
Z(9) w p /B. For V = 50 and p Q = 4.5 x 10~ 2 , S 
is numerically estimated as 4.45 and thus Z(9) ~ 
1.0 x 10~ 2 . For other choice of po, 4.5 x 10~ 3 
and 2.0 x 10~ 3 , po/B becomes 1.0 x 10 -3 and 
4.5 x 10 -4 , respectively. These values agree with 
those of flattening observed in Fig.m and we thus 
find that flattening for these Z(9) is not caused 
by the error in P(Q) but is due to the additional 
term in Eq. (J8J. On the other hand, flattening for 
the po = case in Fig. ^ is due to the errors in 
P(Q) as investigated in (I). 

Having distinguished two kinds of flattening, 
we test whether the MEM can properly repro- 
duce the true one. We employ the case for 
po = 4.5 x 10~ 2 , as an example. Figure displays 
Z^ a \9) for various values of a, which is calculated 
by Eq. (JSJ ■ The default model is chosen to be the 
constant one, m(9) = 1.0. We find that Z^ a \9) 
are independent of the values of a. They show 
clear flattening, which are in agreement with the 
result of the Fourier transform and also with that 
of the Poisson sum formula, Eq. @ ■ The proba- 
bility prob(Z n \P(Q),I,m) shows a narrow peak 
located at a w 15. Thus Z{9) agrees with Z {a \9) 
in Fig. and the MEM yields true flattening. 

With other choice of the value of po, we re- 
peat the same procedure. For p n = 4.5 x 10~ 3 , 
the MEM yields true flattening as in the case 



po = 4.5 x 10~ 2 . This is plotted in Fig. For 
Po = 1.0 x 10~ 3 , however, this is not the case. We 
find that the critical value of po above which flat- 
tening is properly reproduced is po ~ 2.0 x 10 -3 . 
This critical value of po is associated with the 
magnitude of the error in P{Q). 

In summary, we have applied the MEM to mock 
data with true flattening and found that flatten- 
ing is well reproduced. This is in contrast to the 
case of fictitious flattening studied in (I), where 
the MEM could predict no fictitious flattening 4. . 
Applicability of the MEM associated with the 
magnitude of the errors in P(Q) will be reported 
elsewhere. 
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